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I. INTRODUCTION 



Inflationary models predict a nearly scale invariant spectrum of both scalar and tensor 
perturbations, both of which arise from the quantum fluctuations of nearly free fields and are 
Gaussian in character. The scalar perturbations arise from the the quantum fluctuations of 
an inflaton field [TMH(] , and have apparently been observed in the temperature fluctuations of 
the cosmic microwave background |6( . The tensor perturbations arise from the fluctuations of 
quantized linear perturbations of de Sitter spacetime 0-0] , but have not yet been observed. 
In both cases, the Gaussian nature of the fluctuations and the approximate scale invariance 
arise from the properties of free quantum fields. 

Coupling of the inflaton or graviton fields to other fields can modify these conclusions. 
For example, the coupling of graviton modes to the expectation value of the quantum stress 



tensor of a conformal field was recently treated in Ref. [10]. It was shown that graviton 
modes can acquire a one-loop correction which increases their amplitude in a way which 
depends upon the duration of inflation and upon the wavenumber of the mode. This effect 
will tend to lead to a blue tilt to the spectrum of tensor perturbations, but will not change 
their Gaussian character at the one-loop level. 

However, an additional source of perturbations is quantum fluctuations of the stress ten- 
sor. The effects of stress tensor fluctuations in generating density perturbations have recently 



been studied in Refs. [11|, [12J, where a non-Gaussian, non-scale invariant contribution was 
found. Furthermore, this contribution can also depend upon the duration of inflation and 
potentially be used to place limits on this duration. The effect studied in Refs. [ll|, ll2j arises 



from the quantum fluctuations of the comoving energy density of a conformal field in its 
vacuum state. The resulting density perturbations are a non-Gaussian, non-scale invariant 
component to be added to the effect of inflaton field fluctuations [l|-[5j]. 

Fluctuations of other components of the stress tensor are capable of creating tensor 
perturbations. The purpose of the present paper is to address the creation of gravity wave 
fluctuations by stress tensor fluctuations of a conformal field in its vacuum state. These 
can be called passive fluctuations of gravity, as opposed to the active fluctuations discussed 
in Refs. (^[sl- The radiation of gravity waves by stress tensor fluctuations of matter fields 
in thermal states in flat spacetime was discussed in Ref. [l3|. Matter fields in the vacuum 
state in flat spacetime cannot radiate due to energy conservation, but in a time-dependent 
spacetime, such radiation is possible. 

Unless otherwise noted, units in which G — c — H — 1 will be used, where G is Newton's 
constant. 



II. GRAVITATIONAL RADIATION IN AN EXPANDING UNIVERSE 

Here we review the formalism needed to compute gravitational radiation by a time- 
dependent source. We consider a spatially flat Robert son- Walker universe, for which the 
metric may be written as 

ds 2 = -dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ) = a 2 (rf) {-drf + dx 2 + dy 2 + dz 2 ) . (1) 

Here t is the comoving time, and 77 the conformal time. Let 7 MV be this background metric, 
and h^ v be a linear perturbation, 

9p, v = 7mk + hftv • (2) 
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Here we are concerned with tensor perturbations, and impose the transverse, tracefree gauge 
in which 

h^. u = 0, h* = h = 0, and hTuy = , (3) 

where u u = 8" is the four velocity of the comoving observers, and the semicolon denotes the 
covariant derivative on the background spacetime. These conditions remove all of the gauge 
freedom, and leave the two degrees of freedom associated with the polarizations of a gravity 
wave. 

Lifshitz 14| showed that, in the absence of a source, the mixed components satisfy 
the scalar wave equation, 

□V = > (4) 

where 

□ = -4^(^=77^) (5) 

is the scalar wave operator for the metric of Eq. (JTJ). A consequence of this result is that 
gravitons in the spatially flat Robertson- Walker spacetime behave as a pair of massless, 
minimally coupled quantum scalar fields (l5| . 

In the presence of a source, the metric perturbation satisfies an inhomogeneous equation 

□V = -16vrV> (6) 

where S"(x) is the transverse, tracefree part of the stress tensor of the source. It satisfies the 
conditions in Eq. fl3]), and is most conveniently defined in momentum space. The solutions 
of Eq. (j3J) in the spatially flat Robertson- Walker spacetime may be taken to be plane waves 
of the form 

h;(x) = e ; f k ( v )j k -*, (?) 

where fk(v) is a solution of 

+ ,V/ = 0, (8) 
and = e^(k, A) is a polarization tensor which satisfies 

ef = e ; Uu = e ;k u = o. (9) 

If we take vector k to be in the z-direction, then the independent linear polarization tensors 
can be taken to have the nonzero components 

ef = -e/ = ^ , (10) 

for the + polarization, or 

e y = e x = — (11) 
y a/2' 

for the x polarization. 

Denote the spatial Fourier transform of any field A(rj, x) by 

4(17, k) = ^3 jd 3 xe^A( V , X ). (12) 
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In momentum space, the transverse, tracefree part of a stress tensor is defined by the pro- 
jection 

S, U (V, k) = £ e/(k, A) e a "(k, A) T>, k) . (13) 

A 

Thus given a stress tensor T a a(r),^) in coordinate space, we first take a Fourier transform 

to find T a JT), k), then find S^^k) using Eq. ffTB"]) . and finally take an inverse Fourier 
transform to find S"(r), x). 

Let Gr(x, x') be the retarded Green's function for the scalar wave operator, which satisfies 

nG R (x,x') = - S{x ^ ) , (14) 

and Gr(x,x') — if t < t'. Here □ is understood to act at the point x. The gravity wave 
radiated by the source can be written as an integral over the past lightcone of the point 

x as 

hj(x) = 16tt J d 4 x'^/-j(x') G R {x, x') S^{x') . (15) 

The source here could represent either classical matter or quantum fields. In the latter case, 
the average effect of a quantum stress tensor can be described by the semiclassical theory, in 
which the renormalized expectation value (T^) is used as a source in the Einstein equation. 

The effects of a conformal quantum field upon graviton modes in de Sitter spacetime 
has recently been treated in the context of the semiclassical theory (lQ| . It was found that 
there is a correction to the graviton modes which grows with increasing duration of inflation, 



analogous to the effects found in Refs. [Ill, [12[ and to be discussed in this paper. However, 



the effect found in Ref. |l0( comes only from the the expectation value of the stress tensor, 
not from stress tensor fluctuations. If h" is a classical solution of Eq. 01]), then there is a 
correction term h' " . In the case that the conformal field is the electromagnetic field, the 
fractional correction is 

h! v 1 

"°\HSk, (16) 



h v 



57T 



-p ■ 



where l v is the Planck length, H is the Hubble parameter during inflation, and 5* is the 
expansion factor during inflation. This effect grows with increasing S and k, but its total 
magnitude is limited by the requirement that T < 1 for the one-loop approximation to hold. 



III. A FLUCTUATING SOURCE 

Now we consider the case where the source S^(x) is undergoing fluctuations, leading to 
a fluctuating tensor perturbation, h^(x). The correlation function for the perturbation is 

k; ;{x }X ') = (h;(x) W)> - (V(*)XW)} > ( 17 ) 

and that for the source is 

c; /(*,*') = (s;(x) s;{x')) - (s;(x))(s;(x')) . (is) 
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Their relation follows from Eq. (TT5j) : 



K a P a { x , x ') = (16tt) 2 / d 4 X! a/-7(xi) rf 4 x 2 y/-n/(x 2 ) G R (x, x x ) G R (x', x 2 ) C u f(x 1 ,x 2 ) ■ 



(19) 

The spatial Fourier transform of this equation may be expressed as 



K; p a (v, rf, k) = 64(2tt) 8 J d m d m a\ Vl ) a\ V2 ) G( V , m , k) G(rf, m, k) C; p °( Vl , r, 2 , k) , 

( 2 °) 

where C" f(T]i, r/ 2 , k) and G(r], rj', k) are the Fourier transforms of p a (xi,x 2 ) and of the 
retarded Green's function Gr(x,x'), respectively. 

If k is in the z-direction, then the nonzero components of C " p CT (?7i,^2, k) for the + po- 
larization are 

C_ s~i x x y y _ /~i x y _ /~i y x /oi \ 

+ — u x x — U y y — ~ K ~ J x y ~ ~ K - J y x • \ ZL ) 

Similarly, the nonzero components for the x polarization are 

n — n y y — A x x _ A x y _ A y x 

x — ^y y — U y x ~ ^ x y • 

In fact, the stress tensor correlation functions for both polarizations are equal in our case, 
so we may drop the polarization label and write 

C(r] h r] 2 ,k) = C+(r]i,r] 2 ,k) = C x (r] ll r] 2} k) . (23) 

Furthermore, the correlation function C(rji,rj 2 , k) for the conformal field in Robertson- 
Walker spacetime may be related to the corresponding correlation function for the conformal 
field in flat spacetime by a conformal transformation. First consider a classical stress ten- 
sor T " in Robertson- Walker spacetime which is conformally related to T„ v in Minkowski 

spacetime. The spatial components of these tensors are related by T i J = a _4 7~ 3 . The same 
conformal transformation applies to the quantum stress tensor correlation function. Al- 
though the conformal anomaly in the expectation value of a quantum stress tensor operator 
breaks the conformal symmetry, the conformal anomaly for free fields is a c-number which 
cancels in the correlation function. Consequently, we can write 

C{Vi>V2, k) = GT 4 (r?i) cT 4 (772) C M (Vi ~ V2, k) , (24) 

where Cm{vi~V2, k) is the Fourier transform of the Minkowski spacetime correlation function 
for a fixed component of 7~ 3 . As is shown in Appendix [A] , it may be expressed as 

k 5 f 1 

G M (Vi ~ V2, k) = ~ 5127r 5 J du {I - u 2 ) 2 coslkuirjx - rj 2 )} . (25) 
This result applies for either polarization. 



IV. THE POWER SPECTRUM IN INFLATIONARY COSMOLOGY 



The well-known Wiener-Khinchin 16, 17 theorem states that the Fourier transform of 



a correlation function is a power spectrum. A corollary of this theorem is that the power 
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spectrum can normally be written as the expectation value of a squared quantity, and hence 
must be positive. However, the latter result can fail in quantum field theory, and negative 
power spectra are possible. This has recently been discussed in Ref. [HI]. Let C(t — t', x — x') 
be a flat spacetime correlation function. We define the associated power spectrum by a 
spatial Fourier transform at t = t'\ 

P{k) = J^f J d3uetk ' U °^ u )- ( 26 ) 

Given a power spectrum, we can find the correlation function in space at equal times as an 
inverse Fourier transform: 

C(0,u) = J d 3 ke- lku P(k). (27) 

(Here we will use C to denote either a generic or a stress tensor correlation function, and K 
to denote metric correlation functions.) 

The power spectrum for the gravity wave fluctuations is just a spatial component of 
K v cr (i],i]', k) in the limit that rf = rj, and is the same for both polarizations. Thus we 
may combine Eqs. (120]) and (j24p to write the power spectrum at conformal time rj = r] r as 

P(k) = 64(2tt) 8 p d Vl d V2 G(Vr, Vi, k) G( Vr , i] 2 , k) C M ( m - m , k) . (28) 

The possible treatments of the lower limits of integration will be discussed below. We should 
note that the quantity which is usually called the power spectrum in cosmology is not P(k), 
but rather 

V(k) = 4nk 3 P(k). (29) 

It is V(k) which is approximately independent of k for the active gravity wave fluctuations. 
The probability distribution for quantum stress tensor fluctuations is a skewed, hence non- 
Gaussian distribution with non-zero odd moments, although the explicit form has only 



been found in two-dimensional spacetime [19|. Consequently, the gravity wave fluctuations 
produced by stress tensor fluctuations will also be non-Gaussian. 
The Green's function G(r), rj', k) satisfies 

(5 + 4 + V>> < 30 > 

as may be verified by taking a spatial Fourier transform of Eq. (JUJ). Here a' = da/drj. Now 
we wish to specialize to the case of de Sitter spacetime, for which the scale factor is 

aM = -± (31) 

with rj < 0. We may set the scale factor to be unity at the end of inflation, rj = rj r , in which 
case rj r = —1/H. Now Eq. ( |30l) becomes 

(3 + 2Had v + G(i], V >, k) = ■ (32) 
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Comparison of this result with Eq. (71) of Ref. [HI reveals that G(rj, rj', k) differs from the 
Green's function defined in the latter reference by a factor of l/[(27r) 3 a 2 {rj')}. Consequently, 
we may use the result of Ref. to write 

H 2 

G(rj, T]', k) = ^ k y 3 {-k(V- i) cos[k(r] - 77')] + (1 + k 2 777/) sin[k(rj - 7/)] } . (33) 

Next we turn to a discussion of some the possible initial conditions which can be imposed 
on solutions of Eq. ([28]) . 



A. Sudden Switching 

Here we impose the initial condition that the metric fluctuations vanish at rj = i] . The 
power spectrum of tensor fluctuations at the end of inflation, rj = r] r = — 1/H, is then given 
by 

r -i/H r -i/H 

P(k) = P s (k) = 64(2vr) 8 / d Vl / drj 2 G( V , Vl , k) G(rj, rj 2 , k) C M (vi ~ V2, k) . (34) 

Jrin J nn 



The integrals in Eq. f[34|) may be evaluated, using for example the algebraic computer pro- 
gram Mathematica. In the limit that A; 1 770 1 ^> 1, the result is approximately 

P s (k) - -f^f(l + k 2 H- 2 ) . (35) 
air k 

There are several remarkable features of this result: its negative sign, its blue tilt, and 
the fact that it grows with increasing 1 770 1 - The possibility of negative power spectra was 
discussed in Ref. [18| , where it was shown that such spectra arise naturally in quantum field 
theory for the fluctuations of quadratic operators, such as quantum stress tensors. Indeed, 
the power spectrum associated with the fluctuations of the transverse, tracefree part of the 
electromagnetic stress tensor is given by the r]i = 772 limit of Eq. ( T25|) . 

6 " ((U) = -9S^ (36) 

which is negative. Negative power spectra are always associated with coordinate space corre- 
lation functions which are singular in the coincidence limit. This is the case for stress tensor 
correlation functions. They are also associated with the opposite correlation versus anti- 
correlation behavior as compared with a positive power spectrum with the same functional 
form. This means that C(r) changes sign of the sign of P{k) changes, so events that were 
correlated become anticorrelated and vice versa. The spectrum is also not scale invariant, 
and tilted toward the blue end of the spectrum because \V s {k)\ grows with increasing k. 

Another feature of Eq. fl35l) is that the power spectrum for the gravity waves grows as 
rj 2 , which means that it is proportional to the square of the scale factor change between the 



initial time and the end of inflation. This is analogous to the results found in Refs. [11|, [12| for 
the power spectrum of density fluctuations produced by quantum stress tensor fluctuations. 
In both cases, the growth of fluctuations can potentially be used to place upper limits on 
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the duration of inflation, as will be discussed in Sec. [V] The net expansion factor during 
inflation is 

S = H\ m \, (37) 

so we may write Eq. ( 135]) as 

P s (k) = -^(l + eH-*). (38) 

37T d k 

The coordinate space correlation function associated with this power spectrum is given by 
Eq. (J2ZI> to be 

a» = " w t 1 -TP?) ■ < 39 > 

This gives the correlation of points at spatial separation r at equal times. Note that it may 
be either positive or negative. 

Although the power spectrum and the associated correlation function grow with increas- 
ing S or energy scale k, the perturbative approach used here require that |i^(r)| <C 1, which 
places a limit on the magnitude of the effect. 

It is informative to compare the results of this subsection with the flat spacetime limit. 
If we set a = 1 in Eq. (I3"0l) and solve for the flat space Green's function, the result is 

G M ( V - r,', k) = — sm[k( V - rf)] . (40) 

If we use this Green's function in Eq. (128]) . the resulting power spectrum becomes 

^k 

P M (k) = -f- 3 , (41) 

where a rapidly oscillating term which depends upon the integration interval has been 
dropped. The associated coordinate space metric correlation function at equal times is 
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K M(r) = ^- A . (42) 

This function simply describes Planck-scale fluctuations, which are presumably unobserv- 
able. The main point is that the fluctuations do not accumulate in flat spacetime due to 
anticorrelations. In a curved spacetime, such a de Sitter space, this is no longer the case, and 
the anticorrelated fluctuations need not cancel. In the calculations, the crucial difference is 
between the flat spce Green's function, Eq. (J4"0|) . and that in de Sitter space, Eq. (133]) . 



B. Exponential Switching 

In the previous subsection, the interaction between the quantum stress tensor and the 
gravitational field was taken to be switched on suddenly at rj = t] Q . One might be concerned 
that either the sign of P(k), or its growth with increasing \t] Q \ are artifacts of this sudden 
switching. Here we investigate a model in which the interaction is switched on gradually. 
We replace the step function 9(rj — t]q) by an exponential function, e pv , with p > 0. This 
function vanishes as rj — > — oo, and in the limit of small p, is close to unity by the end of 
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inflation. The effect of this switching function is effectively to switch on the interaction on 
a conformal time scale of order 1 770 1 5 where 770 = — 1/p. Equation (1341) is replaced by 

r-l/H r-l/H 

P e (k) = 64(2tt) 8 / drix s e (r) X ) I dT] 2 s e (r) 2 ) G(r], 771, k) G(rj, rj 2 , k) C/Zat^i - V2, k) , 
J —00 J —00 

(43) 

where the switching function is 

s e (v) = e m ■ (44) 
In the limit of small p, Eq. fj4*3]) leads to 

m) H .?^m + 0(M = jm , ,45) 

8n z k z p m z k z 

where 5, given by Eq. ( 1371) . is the expansion between t] = t]q = —1/p and the end of inflation. 

Again the power spectrum is negative, blue tilted, and grows with increasing S although 
now linearly. In this case, the equal time spatial correlation function is an inverse Fourier 
transform of P e (k) given by 

K e {r) = -^. (46) 

Here we have dropped a delta-function term proportional to <5(x), which will not contribute 
to measurements made at distinct spatial locations. Note that because a{rj) = l/(H\r}\), if 
the switching time A?? is of order |?yo|, then the scale factor approximately doubles during 
the switch-on. For example, a(rj/2) = 2a(rj). In terms of comoving time t, where a(t) = e Ht , 
this corresponding to a time interval of At w l/H, or one horizon crossing time. Thus the 
switch-on time in this model is of order of the horizon crossing time. 



C. Adjustable Width Switching 

The switching function s e {rf) used in the previous subsection contains only one parameter, 
p, which regulates both the effective duration of inflation and the period over which the 
switching occurs. It is instructive to consider a more general function with two parameters: 

8 ™W= i + Jno-v)/* - (47) 

This function, analogous to the Fermi-Dirac distribution function, changes from zero to 
unity when r\ w r/ over a time scale of Ai] w a. The resulting power spectrum, P aw (k), is 
given by Eq. (J4"3"j) . with s e (rj) replaced by s aw (rj), and may be expressed as 



PaM = "2^3 du (1 - U 2 ) 2 (I 2 C + Is) ■ 



(48) 



Here 

Ic 

and 



/•oo 

/ dx g(x, x r ) s(x) cos(ux) , (49) 

J x r 

POO 

/ dxg{x,x r ) s{x) sin(w:c) , (50) 

J X r 
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with 

g(x, x r ) = (x — x r ) cos(x — x r ) — (1 + xx r ) sin(x — x r ) . (51) 

We use the notation, x = —kr], x r = —krj r , and s(x) = s aw (v)- The dominant contributions 
to the integrals in Iq and Is come from values of x of order xq = —kr] 3> x r , so we may 
write 

g(x, x r ) ~ x(cos x — x r sin x) . (52) 

The resulting integrals may be evaluated using Eqs. (1B8[) and (1B9j) . which are derived in 
Appendix IB1 The final result, when a > 1/k, is 

^„(A:)^--^(1 + A; 2 F- 2 ). (53) 
27r A;" 1 a 

Apart from numerical factors, this result contains both Eqs. ( 135]) and (|45j) as special cases. 
If a f=s 1/fc, then we return to the sudden switching case of Eq. (1351) . On the other hand, 
if a fa |?7o|, we find Eq. (H5|) . up to numerical constants. The fact that the constants do 
not match exactly may be due to the approximation used in deriving Eqs. (1B8|) and (1B9j) 
(g <C 1) not being very good near u = 0. In summary, if Arj < 1/fc, we obtain the sudden 
result, Eq. (13"5|) . proportional to S 12 , and if Ar] ps \r) \, we obtain Eq. fH5|) . proportional to 
S". Intermediate switching times lead to Eq. (IBUj) . 



V. IMPLICATIONS OF THE POWER SPECTRUM 



A. Initial Conditions and the Transplanckian Issue 

Although the results in the previous section depend somewhat on the rate at which 
the coupling between the quantum stress tensor fluctuations and the gravitational field is 
switched on, in all cases the power spectrum grows as a power of S, the expansion from 
the initial time to the end of inflation. Thus we need an interpretation which suggests a 
reasonable value for this time, \tjq\. One possibility is to take this time to be the onset 
of inflation. This imposes the initial condition that the gravity wave perturbations vanish 
at the beginning of inflation. In this case, S becomes the total expansion factor during 
inflation. A possible objection to this interpretation is that it can lead to contributions from 
transplanckian modes. This raises the question of whether our perturbative treatment can be 
trusted, as relations such as Eq. f|T5l) are lowest order approximations in the dimensionless 
coupling constant (£ p k) 2 . The transplanckian issue has been extensively discussed in the 
contexts of the Hawking effect and of cosmology. Hawking's original derivation [20] of black 



hole radiance relies upon modes which begin far above the Planck energy. The fact that the 
Hawking effect gives a beautiful unification of gravity, thermodynamics, and quantum theory 
can be considered to be a powerful argument to take transplanckian modes seriously. It is 



true that it is possible to derive the Hawking effect without transplanckian modes [21|, [22 
but only at the price of introducing modified dispersion relations which break local Lorentz 
symmetry and hence postulate new physics. There has been an extensive discussion of 
the possible role of transplanckian modes in inflationary cosmology. (See Ref. [l2| for a 
lengthy list of references.) The effect discussed in this paper has the potential to serve as 
an observational probe of transplanckian physics. 
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There is an alternative possibility [l2j], which is to take the initial time at which the 
perturbation vanishes to depend upon the mode, and to be the time at which a given mode 
redshifts below the Planck scale in the comoving frame. This avoids the transplanckian 
issue, but at the price of introducing a non-local and frame dependent prescription, which 
is analogous to introducing non-Lorentz invariant modified dispersion relations. In the 
remainder of this paper, we will explore the consequences of adopting the former prescription 
whereby S is the total expansion factor during inflation. 

The dependence of the gravity wave spectrum upon a positive power of S might seems to 
contradict a theorem due to Weinberg 23[, which was generalized by Chaicherdsukal [24 



This theorem states that radiative corrections during inflation should not grow faster than a 
logarithm of the scale factor. However, as was discussed in more detail in Ref. [l2l ]. density 
perturbations which are proportional to a power of S are really due to high frequency 
modes at the initial time, and are hence always large rather than growing. This comment 
also applies to the effects found in Ref. [Toj and in the present paper. 



B. Numerical Estimates 

We may use the coordinate space correlation functions, K s (r) and K e {r) , to estimate the 
physical effects of the gravity wave fluctuations on various scales. However, these functions 
describe the primordial fluctuations at the end of inflation. After the end of inflation, modes 
which are outside the horizon remain approximately constant until they re-enter the horizon. 
(For a more detailed discussion, see, for example, Ref. (25[.) After that point, they redshift 
with their amplitude proportional to I /a. Let a# c be the value of the scale factor at which 
a mode associated with coordinate length r re-enters the horizon, and a now be the present 
value of the scale factor. The present value of the correlation function is then 

K now (r) = K(r) (^-) . (54) 

For the sudden switch model, this becomes 

4H*S*e>„ / 2 \ (a Hc \> 



and for the exponential switch model it is 



H 3 S> 



K now Mr) = --^p 1 ( — ) • (56) 

Here the factors of the Planck length i v are written explicitly. 

Let Er be the reheating energy at the end of inflation. This energy has since been 
redshifted to that of the cosmic microwave background. We set a = 1 at the end of inflation, 
so that 

anow ~ 2.5 x 10- W ' 1 ' 

The proper length scale today associated with coordinate distance r is 

t = a n ow r . (58) 
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We assume that reheating is efficient, so the vacuum energy at the end of inflation is of order 
E R , and 

H 2 = ff p E R . (59) 

We also assume that the scale of interest was outside the horizon at the end of inflation, so 
that Hr > 1. We may combine all of these results to write 

|A_ S | = 10- ( %V f TtJ^V OY S 2 , (60) 



and 



10 16 GW 7 V a, 



l^-l = (»] If^l 5- (61) 
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Let us first consider the case of perturbations of the order of the present horizon size, 
t m 10 61 i p . In this case, a# c ~ a now . Data from the WMAP satellite constrain these 
perturbations to satisfy h < 10~ 5 , so that |if n£W | < 1CT 10 . Consequently, the sudden switch 



model leads to 



10 16 GelA 



3 



5< 10 34 (^-^j , (62) 



and the exponential switch model to 



W 16 GeV\ 
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S<10 W ^— —j . (63) 

These constraints on the total expansion during inflation are compatible with adequate 
inflation to solve the horizon and flatness problems, S > 10 23 . Because K < 0, quantum 
stress tensor fluctuations during inflation will tend to produce anti-correlated gravity wave 
fluctuations. Note that in this example, |i^ nom | ~ K{r) < 1CT 10 , so the criterion for the 
validity of the perturbative calculation, |iT(r)| 1, is satisfied. 

Now we wish to consider perturbations which are well within the present horizon. For 
this purpose, we need an approximate model for the current matter content of the Universe. 
Although the dominant component today is the dark energy, this is likely to be a recent 
phenomenon. If the dark energy is due to a cosmological constant term, it does not redshift 
and hence does not grow as we go backwards in time. Here we assume that the Universe 
was radiation dominated, a oc t 1 / 2 , for t < t eq and subsequently matter dominated, a oc 
Furthermore, we assume 



a 



10" 4 , (64) 

^now 

so that teqjtnow ~ 1CU 6 . A perturbation with proper length t enters the horizon at t = £# c = 
£. If we assume that £ < t eq , then we may write 

a ^ 2 ~l(T 63 ^. (65) 



f 



If we insert this relation into Eq. f l5D|) . the result is 



|A ,_,^ M M (^)V m 
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In the case of the exponential switch model, the factors of i cancel, 

^-^^(i^)' 5 ' < 67) 

leading to a scale independent correlation function on scales i < 10 23 cm. 

If the magnitude of these fluctuations is sufficiently large, they should produce background 
noise in gravitational wave detectors, which has not been observed. LIGO has placed lim- 
its 26|] of h < 10 -24 on scales of the order of 10 2 km, corresponding to |-fT ncm ,| < 10~ 48 , and 
leading to the constraints 

10 10 GeV\ 3 



S < 10- (— — j . (68) 
for the sudden switch model, and 

s< (69) 



for the exponential switch model. However, these results are compatible with adequate 
inflation to solve the horizon and flatness problems only if 

E R <10 10 GeV (70) 

for the sudden switch model, and 

E R <10 13 GeV (71) 

for the exponential switch model. In this example, |-K"(r)| ~ 10 23 K now < 10 -25 , so again the 
requirement that \K\ be small is fulfilled. 



VI. SUMMARY AND DISCUSSION 

We have seen that quantum stress tensor fluctuations are capable of creating gravity waves 
during inflation. The resulting spectrum has several properties, including negative power 
and an amplitude which grows with increasing duration of the inflationary period. Negative 



power spectra, although forbidden by the Wiener- Khinchin theorem 16j, ll7|, can arise in 
quantum field theory [18| , especially in quantum stress tensor fluctuations. A negative power 
spectrum can be viewed as interchanging correlations and anti-correlations, as compared to 
a positive power spectrum of the same functional form. 

We find that the amplitude of the gravity wave spectrum is proportional to a positive 
power of S, the change in scale factor during inflation. A similar dependence was also found 
in Refs. [ill. 12| . for the effects of stress tensor fluctuations on density perturbation and 
in Ref. [l0| for the correction to gravity wave modes from expectation value of the stress 
tensor of a conformal field. The gravity wave power depends somewhat upon the details of 
the initial conditions, being S 2 if one integrates the equations directly from a state of zero 
fluctuations, and being S if the interaction between the fluctuating matter stress tensor is 
supposed to be switched on over a finite interval of the order of the horizon size in comoving 
time. In all cases, the primordial power spectrum of gravity wave fluctuations is negative 
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and greater in magnitude at shorter wavelengths. This non-scale invariant spectrum of 
fluctuations will be highly non-Gaussian, due to the non-Gaussian character of quantum 
stress tensor fluctuations. 

These gravity wave fluctuations are potentially observable. Longer wavelengths could 
alter the polarization of the CMB, and be detected in the same way as the active fluctuations. 
Shorter wavelengths could potentially be detected by Earth or space based gravity wave 
detectors. The fact that they have not yet been detected can be used to infer constraints on 
the duration and energy scale of inflation. 
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Appendix A: Flat Space Stress Tensor Correlation Functions 

In this appendix, we derive the explicit expressions for the flat space correlation functions 
utilized in Sect. 1111} especially Eq. (125]) . All of the expressions in this appendix refer to flat 
spacetime, so here we drop the subscript "M". We may use the results of Ref. 27|, where 
the electromagnetic field stress tensor correlation function was shown to be 

C^\x,x') = 4(d^D)(d a d x D) + 2g^(d a d a D)(d x d a D) + 2 g aX (d,d a D) (d„d a D) 
~ 2g^(d u d a D)(d x d a D) -2g va (d^d a D)(d x d a D) 
~ 2g uX (d fl d a D)(d a d a D) -2 9llX (d u d a D)(d a d a D) 

+ {g^QuX + 9ua9ti\ - g^9ax) (d P d a D) (d p d a D) . (Al) 

Here 

D = D(x-x')= 1 (A2) 

4:7T Z {X — X ) Z 

is the Hadamard (symmetric two-point) function for the massless scalar field. For our 
purposes, it is sufficient to compute a single component, such as C xyxy . The result is 

3 

^ = 7r 2 [(t _ t0 2_ r 2 ] 4 > ( A3 ) 

where r = |x — x'| and r = t — t' . The spatial Fourier transform of this expression is 
or equivalently, 

C*v*v( T , k) = r / du{l- u 2 ) 2 cos(kur) , (A5) 



5127T 5 

14 







which is Eq. (|25|) . This expression may be verified by checking that 

C xyxy^ r ) = J d 3 k e - ik '( x - x ') C™(r, k) . (A6) 

Appendix B: Evaluation of Sampling Function Integrals 

In this appendix, we will evaluate some of the integrals needed in Sect. II VQ which 



involve the function s aw , defined in Eq. (H7|) . We begin with expression 3.411.23 in Ref. [28 
which states that 

dx = — 7T 2 CSc(7T/i) COt(7T/i) , (Bl) 



-oo 1 + 

for < Re(fi) < 1. This implies that 



l + e 



- dx = IT csc(7r^t) . (B2) 



This may be verified by taking a derivative of Eq. (1B2p with respect to fi, and by noting 
that when /i = 1/2, this relation becomes 

oo x/2 roo ^ 

dx = 2 - dy = 7T . (B3) 



-oo 1 + e x Jo 1 + y 

This confirms that there is no additional constant in Eq. (IB2p . Next we may take the limit 
in which fi — > iq to write 

dx = . (B4) 

_j 1 + e x e 7 ^ - e _7r 9 

However, we need integrals over a semi-infinite range of the form 



oo piqx 

dx = be iqX0 



I _|_ e (x-x )/b 



oo ^tqbz r — i — e iqbz 



-no 1 + e z / ^ l + e 



dz — dz 



(B5) 



The second integral on the right-hand side of the above equation may be approximated by 
setting the denominator of the integrand to unity: 



— 5 — p iqbz ■ 

dz « — e iq{ " o) + O(e- xo/b ) . (B6) 

-oo 1 + ^ qb 

Thus, 

iqx n ■ iqx ■ 

dxP z-J?ll + - e ^. (B7) 



X -)- e (x-x )/b e nq _ e -irq g 

If we take a derivative with respect to q, then the real and imaginary parts of the resulting 
expression become, for qb ^> 1, 

f°° x sin ox , , . . 

J 1 + e (»-x )/6 dX W " 27rX ° 6 C ° S ^ X °) 6 ( B8 ) 
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and 



f x cosqx , . . _ nh , 

/ -, r-rr- ax « 27tx o smfoxo) e y . B9 

y a 1 + e( x - x °y b K J 
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